Spin-orbit effects in nanowire-based wurtzite semiconductor quantum dots 
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We study the effect of the Dresselhaus spin-orbit interaction on the electronic states and spin 
relaxation rates of cylindrical quantum dots defined on quantum wires having wurtzite lattice struc- 
ture. The linear and cubic contributions of the bulk Dresselhaus spin-orbit coupling are taken into 
account, along with the influence of a weak external magnetic field. The previously found analytic 
solution for the electronic states of cylindrical quantum dots with zincblende lattice structures with 
Rashba interaction is extended to the case of quantum dots with wurtzite lattices. For the electronic 
states in InAs dots, we determine the spin texture and the effective g-factor, which shows a scaling 
collapse when plotted as a function of an effective renormalized dot-size dependent spin-orbit cou- 
pling strength. The acoustic-phonon-induced spin relaxation rate is calculated and the transverse 
piezoelectric potential is shown to be the dominant one. 

PACS numbers: 73. 21. La, 71.70.Ej, 73.61.Ey, 72.25.Rb, 



I. INTRODUCTION 

Semiconductor nanowires catalytically grown from 
nanoparticles, also called nanorods and nanowhiskers, 
provide a promising platform for spintronic devices. One 
of their advantages is that small quantum dots can be 
conveniently defined in these nanowires by employing dif- 
ferent techniques. For example, quantum dots can be 
obtained by varying the composition during the growth 
process^^ in the same way that quasi-two-dimensional 
superlattices are grown using molecular-beam epitaxy 
or metalorganic chemical vapour deposition techniques. 
This fabrication method produces sharply defined quan- 
tum dots with square-well confinement in the longitudi- 
nal direction with highly controllable lengths. This type 
of quantum dot has been thoroughly studied in the last 
ten years and its spin-related properties have attracted 
much interest. In particular, the effective g-factor in InAs 
nanowire-based quantum dots has been measurecP for 
the first few electrons entering the dot and a strong de- 
pendence on the dot size has been established. More- 
over, studies of spin-orbit interaction and spin relaxation 
have been made in InAs nan owire-based double-quantum 
dots with two electrons!^ A less precise but very flex- 
ible type of quantum-dot structure can be obtained by 
electrostatic confinement given by thin metallic gates de- 
posited perpendicularly to the wires, typically made of 
InAs or InSbP^ The gate voltages, which are gener- 
ally time-dependent, allow for the definition of the quan- 
tum dots as well as for the control of charge transport 

and spi n manip ulation via electric-dipole spin resonance 
(EDSR)P™ 

An important aspect of semiconductor nanorods, and 
also of the quantum dots defined in them, is that they 
often display the wurtzite crystal structure even though 
the constituting material has a zincblende structure in 
the bulkP^ This structural change has important conse- 
quences for the spin properties of the dots. Notably, while 
in the zincblende semiconductors the leading Dresselhaus 



spin-orbit coupling term is cubic in k, in wurtzite crystals 
a linear term appears.^ A realistic theoretical study of 
the spin-orbit effects in wurtzite quantum dots has now 
become possible due to recent advances in the character- 
ization of the band-structure of different wurtzite mate- 
rials. While the linear-in-k Dresselhaus contribution to 
the conduction band energy of wurtzite crystals has been 
known for some timej^the cubic term in k has been ob- 
tained only recently within the k • p approximation for 
different wurtzite semiconductors {HEU pje and Pryorpf 
on the other hand, calculated the band-structure param- 
eters of several binary compounds which normally dis- 
play the zincblende structure in the bulk assuming that 
they have the wurtzite structure, in order to make them 
available for studies of wurtzite nanowires made of those 
same materials. These new data pave the way to a realis- 
tic study of nanowire-based quantum dots with wurtzite 
structure, which we undertake here. 

Disk-shaped quantum dots with the Rashba structural 
spin-orbit interaction have been extensively studied the- 
oretically and have been shown to admit an analyti- 
cal solution for their energy eigenstates, without 1 ' and 
with^ an applied perpendicular magnetic field. This so- 
lution, as we show here, can be conveniently extended 
to wurtzite quantum dots having cylindrical symmetry 
around the crystal c-axis, either with flat ( "disk" ) or elon- 
gated ( "rod" ) geometry. This is the case since the linear 
wurtzite Dresselhaus coupling is mathematically equiva- 
lent to the Rashba linear spin-orbit coupling character- 
istic of asymmetric semiconductor quantum wells, and 
furthermore, the newly obtained cubic term of wurtzite 
admits in the quasi-two-dimensional case the same eigen- 
states as the linear term. In this work we exploit these 
similarities in order to give solutions of the eigenvalue 
problem of the wurtzite quasi-two-dimensional structures 
and cylindrical quantum dots. In a confined geome- 
try, the wurtzite cubic term of the Dresselhaus coupling 
gives rise to an additional linear contribution that rein- 
forces or counteracts the bare linear term. As we will 
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see below, this reinforcement can be actually much big- 
ger than the original linear term, opening up an unex- 
plored regime of strong "Rashba-like" spin-orbit coupling 
in quantum wells and dots. Also, new possibilities appear 
when these linear Dresselhaus terms are combined with 
the standard Rashba term due to structural asymme- 
try. Indeed, further reinforcement or cancellations could 
possibly be achieved by tuning the symmetry and di- 
mensions of the structure with the help of an external 
gate voltage. Thus, flexible schemes of spin-orbit cou- 
pling cancellation could be implemented leading to very 
long spin relaxation times in wurtzite structures having 
particular geometric shapesP^The Dresselhaus term, en- 
abled by the bulk inversion asymmetry, has been shown 
to yield the dominant coupling mechanism in cases of 
important structural asymmetry, like that of an extrinsic 
impurities giving rise to the impurity band of n-doped 
semiconductors p3 

The article is organized as follows. In Section [TT] we in- 
troduce the effective Hamiltonian including the spin-orbit 



coupling for wurtzite structures. In Section III we obtain 
the electronic states of quasi-two-dimensional structures 
considering first only the linear spin-orbit coupling (Sub- 
section III A I and then the full Hamiltonian with the cu- 



bic spin-orbit term and the Zeeman energy (Subsection 



III B I . In Section IV the solution of the previous Section 



is used to solve the problem of thin cylindrical quantum 
dots with hard-wall confinement potential. In Subsec- 
tion |IV A| we present the general analytical solution of 
this problem and the energy levels calculated numeri- 
cally. In Subsection |IVB| we explore the spin structure 
of the one-particle eigenstates. The experimentally ac- 
cessible effective g-factor of the quantum dots is studied 
in Subsection |iV C[ and in Subsection |IV D| wc discuss the 
spin relaxation due to the coupling to phonons. Section 
[V] provides concluding remarks. 



II. INTRINSIC SPIN-ORBIT COUPLING IN 
WURTZITE-BASED CONFINED GEOMETRIES 

Within the envelope-function approximation for 
conduction-band electrons in wurtzite semi conductors, 
the effective quantum-d ot Ha miltonian^ 2 ^ incorporat- 
ing the linear^ and cubic^^ Dresselhaus spin-orbit cou- 
plings reads 



H — H + Hi + H 3 + H z , 



H = ~ — r + V c (x,y,z), 
2 m* 



Hi = a (k y o x - k x a y 



(1) 

(2) 
(3) 



where V c is a nanoscale confinement potential, cr is the 
spin operator, a, 7, and b are material-dependent pa- 
rameters, m* is the effective mass and g* the bulk ef- 
fective gyromagnetic factor, /iB is the Bohr magneton, 
and B is an external magnetic field assumed to be ap- 
plied in the z-direction that coincides with the c-axis 
of the wurtzite structure. Here we include the mag- 
netic field only through a Zeeman term since we will 
consider only relatively weak fields whose orbital effects 
can be safely ignored. In what follows we will con- 
sider quasi-two-dimensional and cylindrical quantum-dot 
structures. Catalytically grown nanorods made out of 
materials which have the zincblende crystal structure in 
the bulk can adopt either the zincblende or the wurtzite 
structure depending on the size of the nanoparticle seed 
and other growth conditions. Experimental data allow- 
ing to determine a, 7, and b are not yet available, so in 
our study we will rely on the theoretical estimates ob- 
tained by De and PryorP^l Motivated by anticipated ap- 
plications to nanowires, these authors calculated all the 
relevant band-structure parameters assuming a wurtzite 
structure for the semiconductor binary compounds that 
have a zincblende structure in the bulk. An asymmetry 
in the z-confinement would add a Rashba term, resulting 
in a renormalization of a. In order to give a wider appli- 
cability to our results, whenever possible we will present 
them for reasonably large ranges of parameters so that 
they can be adapted to different materials and to param- 
eters newly obtained, experimentally or theoretically. 



III. QUASI-TWO-DIMENSIONAL SYSTEMS 

Before tackling the quantum-dot problem it is use- 
ful to consider the eigenvalue problem of a quasi-two- 
dimensional system. Thus, we choose V c — V c (z) which 
confines the electrons only along the z-direction, such 
that H can be separated as H — H^ y + Hq, with an 
in-plane term H^ y = (p x + Py)/2m* and a longitudinal 
part defined by H% — p 2 z /2m* + V c (z). 



A. Linear term 

If we leave aside for the moment the cubic term H3 and 
the external magnetic field, we are left with a situation 
mathematically analogous to the classic Rashba problem 
in which the spin-orbit coupling originates from an asym- 
metric extrinsic potential. Since the total Hamiltonian is 
separable we can start working with the two-dimensional 
problem in the (x, y) plane given by H 2 d = H^ y + H\. 
Its well-known solution is^ 



H 3 = 7 (bk 2 z -k% - k 2 y ) {k y a x - k x a y ) , (4) 



Ck,(r) 



1 



-i [<fk- 

1 



(6) 



Hz = -g*n s Ba z , 



(5) 



E(k,s) 



2 1,2 



K Z k 

2m* 



sak. 



(7) 
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In these expressions and henceforth, A is the area of the 



sample, r = (x,y), k = (k x ,k y ), k = yjk x + k 2 , and ip k 

is the angle of k in polar coordinates. The spin quantum 
number s = ±1 denotes spin-up and spin-down eigen- 
states with respect to the spin quantization axis which 
lies in the xy-plane and is perpendicular to k with a po- 
lar angle ipk — 7r/2. Note that the spin-orbit term in the 
energy has a minus sign compared to the usual Rashba 
expression, coming from the minus sign used in Eq. Q. 
The states (|6| are degenerate for given k and s. This 
plane-wave solution is convenient in most contexts and 
has the advantage that its spin quantization direction 
is position independent. However, ^ does not profit 
from the fact that the z-component of the total angu- 
lar momentum J z commutes with the Hamiltonian and 
therefore provides a good quantum number, which is an 
extremely useful property when one tackles cylindrically 
symmetric nanostructures. The common eigenstates of 
H 2 d and J z are given bjff^ 



Xm,k,s(r,(fi) = 



J m (kr) e lm * 
sj rn+1 (kr) e^ m+1 ^ 



(8) 



The states (|Sj) are degenerate with those of ^ for given 
k and s, and can be expressed as superpositions of them. 
Note that while the spin of the basis states (6j) lies always 
in the xy-plane, that is not the case for the states pi), 
which are superpositions of the states Q within degen- 
erate subspaces. Furthermore, the spin direction in the 
latter is space-dependent while in the former it is not. 



B. Cubic term 

Let us now include the cubic-in-fc term of the Hamil- 
tonian, i?3, given in Eq. Q, and the Zeeman energy, 
Eq. IH]). As usual, we work in the envelope-function ap- 
proximation where the Hamiltonian H is expressed by 
replacing k by — iV. We adopt cylindrical coordinates 
(r, (p, z) and for the in-plane coordinates we have 



H = — 



2m* 



dz 2 



b - 



dz 2 



+ V c (z) + H 1 



Hi + H z . 



(9) 



where the symbol V 2 is used to represent the two- 
dimensional Laplacian. 

Assuming that V c (z) is an infinite potential well of 
length L, the proposed solution of the Schrodinger equa- 
tion H£ = E£ is 



(10) 



£ nTO (r, <p, z) = ip nm (r, <p)Jj- sin (J^j^j 



Uy. 



ip nm (r, f) 
,(r) and 



( r \ pinup 
my ) e 

(r) e i ( m + 1 V 



(11) 



(r) are real functions and 



where 

"0nm ( r i v) is an eigenstate of J z with eigenvalue j z 
m +1/2. The corresponding total energy is 



E„, — E„ 



El 



(12) 



with the radial part E n and the longitudinal energy E* = 
(h 2 /2m*)(nir / L) 2 that is due to the confinement in the 



z-dircction. Plugging ( 10 ) into the Schrodinger equation 
we obtain for u nrn and v nm the equations 



/ m + 1 d \ 

{-^In + h ) U nm (p) + (a' n + YV m ) ( h ^ ] V nm (p) = £ n U nm (p) 

— - -Q - \ U nm (p) = EnVnm(p) 



(13a) 

(13b) 



where 



id 

pdp 



dp 2 



If 



(14) 



In Eqs. ( 13 1 we have introduced R, a parameter to be de- 
fined in the quantum dot context, and u E — h 2 /2m* R 2 , 
as units of length and energy, respectively. This al- 
lows us to define the dimensionless parameters p = r/R, 
K = kR, 7' = 7/-u E J? 3 , and h = gpBB/2u E . The depen- 
dence on the ("longitudinal") quantum number n has 
been incorporated to the in-plane problem via the redef- 



r 



inition of the coupling constant a written in the dimen- 
sionless form 



a + 76 ( ■ 



^7ry 
L~) 



/u e R, 



(15) 



and the in-plane dimensionless energy is given by e n 
E n /u E . 



To solve Eqs. ( 13 ) we make the ansatz 



Unm(p) = J m (Kp), v nm (p) = d n J m+ i(Kp). (16) 
Using the well-known properties of the Bessel functions^ 
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m d 

P dp 

m + 1 d 

P dp 



J m {Kp)=KJ m+1 {Kp), (17a) 
J m+1 (Kp)=KJ rn (Kp), (17b) 



one obtains from ( 13 1 the eigenvalue equation 



K 2 + h - e n a' n K - 7 'K 3 \ / 1 
K a' n K-~f'K 3 K 2 -h-e n ) \d r , 

whose solutions arc 



= 0, 



e n± = K Z ± \/K 2 K-y^r 
Then, the total energy is given by 

(e n ± + e z n )u E , 



h 2 . 



E n± 



(18) 



(19) 



(20) 



with e z n = E^/u E — (nirR/L) 2 , and the corresponding 
wave functions are 



1pKnm{p,<p) 



J m (Kp) e lm * 
d n± J m+1 (Kp) e ^ m+1 ^ 



with 



K 2 



a' n K - j'K 3 



(21) 



(22) 



The obtained solution, Eqs. (19 1 and (21), reduces to 
the one of the linear Hamiltonian analyzed in Sec. |III A[ 
given by Eqs. ^ and when the cubic term and the 
Zeeman energy are neglected. As in the linear Rashba- 
like problem, there are two possible energies e n ± for a 
given value of K. The energies e n ± can be expressed as a 
function of K 2 and are thus independent of the sign of K. 
Because of the (anti-)symmetry of the Bessel functions 
with respect to a change of sign in the argument, the wave 
functions corresponding to ±K are not independent. We 
therefore keep only positive values of K. 

In the presence of a magnetic field, and for in-plane en- 
ergies e n close to zero, Eq. ( 19 1 has solutions with imag- 
inary K = ik. Since J m (inp) = i m I m (Kp), where I m 
is the modified Bessel function of the first kind of order 
m, the corresponding wave functions grow exponentially 
with increasing p and are thus not normalizable in an 
infinitely large system. Such solutions are therefore dis- 
carded in the context of two-dimensional systems, but 
they will become relevant for the case of quantum dots 
discussed in Sec. IIVI 

In Fig. [I] we present (solid lines) the dispersion relation 



( 20 1 for InAs with the parameters suggested in Ref. 
from band-structure calculations (we label them with an 
index r) 

a x = 0.571 eVA 7r = 571.8 eV A 3 , 

b — 4 and an effective mass m* = 0.026 m e Pj The two en- 
ergy branches are plotted for B — T (left) and B = 20 T 
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FIG. 1. Energy dispersion [Eq. |20| ] with (right) and with- 
out magnetic field (left). For subband n = 1, three cases 
are considered: thick solid lines correspond to the full SOC, 
dashed lines to an intermediate case with no cubic-in-k SOC, 



but with the a parameter renormalized by 7 ( 15 1, and dotted 
lines to the bare Rashba-like interaction, linear in k. Blue 
(red) lines correspond to e n - (£n+). The thinner solid lines 
are the lower branches of subbands n = 2 and 3. Inset: Same 
three cases of SOC dispersion relation without the parabolic 
contribution for n = 1. The Zeeman effect and an avoided 
crossing are more clearly distinguished on this energy scale. 



(right). (In this figure we consider a large value of the 
magnetic field with the sole purpose of illustrating more 
clearly its effects on the energy levels.) We also show the 
effect of suppressing the cubic term, but keeping the con- 
tribution of 7 in Eq. ( 15 ) on the linear term (dashed line) , 
as well as the usual Rashba-like case obtained for 7 = 
(dotted line). Blue (red) lines correspond to £„_ (£«+)• 
Thick lines correspond ton = 1, as indicated between the 
two panels. Also shown are the curves of £2- and £3- 



including the linear and cubic spin-orbit contributions 
(thin lines) . For n = 1 and B = T there is a crossing of 
the two branches at K = yV/V (K = 19.34 in our plot). 
An analogous feature present in bulk wurtzite semicon- 
ductors has been discussed in the literature as a possible 
opportunity to implement long-lived spin qubits.^ The 
crossing becomes avoided for finite B, although the level 
splitting can hardly be seen on the right panel of Fig. [T] 
For this reason we plot in the inset the energies subtract- 
ing the trivial parabolic contribution. This allows for a 
smaller energy range such that one can clearly observe 
the Zeeman splitting at K — and the avoided crossing. 

In Fig. [TJ the thin solid lines are the lower branches 
of subbands n = 2 and 3. Even though they lie at suffi- 
ciently high energies so as not to affect our further anal- 
ysis, which concentrates on low energies, we note that 
they could become relevant if the region of the avoided 
crossing mentioned above is explored. Also, we point 
out a potentially interesting degeneracy point of all the 
lower branches of the different subbands, which happens 
at K — l/b-f' (K — 17.58 in our plot), where the curves 
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become independent of n. This massive degeneracy is 
due to the renormalized linear spin-orbit term. Although 
this feature may be physically relevant, we mention that 
higher values of n correspond to higher k z and eventually 
the energies of Eq. ( 20 ) obtained in third-order pertur- 



bation theory in wave-vector cease to be reliable. 



IV. QUANTUM DOTS 

A. Effect of spin-orbit coupling on the energy levels 

We now consider cylindrical quantum-dots with hard- 
wall quantum confinement having radius R and length 
L. The discrete eigenenergies and states of this problem 
will be obtained from the quantum-well solutions found 
in the previous Section. In order to get the energetically 
lowest states, we keep only the lowest subband, n = 1, 
and omit the subindex n from now on. In all cases we 
work with k low enough to stay in the regime of validity 
of the expansion of the effective SOC Hamiltonian up to 



3 rd order in k. 

The eigenstates of the disk-shaped quantum dot have 
to satisfy the circular boundary condition (the hard- 
wall confinement forces a zero of the wave function at 
the dot boundary). This can be achieved at particular 
values of the in-plane energy e for linear combinations 
= c&ipK a m + Cbil>K b m of two degenerate eigenstates 
of the quantum-well problem. Those quantized energies 
are then the eigenenergies of the quantum dot. In the 
general case including a finite magnetic field, there are 
three energy ranges (see Fig. [T]) with different situations: 

i) energies in the low "belly" of the e_ branch, e < — |ft|; 

ii) energies above the energy gap caused at K = by 
the Zeeman splitting, e > \h\; iii) energies in the Zeeman 
gap, — \h\ < e < \h\. We now consider these three cases 
separately. 

Case i) e < —\h\; two real values of K, noted K a and Kb, 
associated to the £_ branch (in the "belly" region) are 
involved in the dot solution. The in-plane wave function 
is thus written as 



j m {K hP y m * 



(23) 



with the hard-wall boundary condition 9 m (p — IjV) = 0. A non-trivial solution (if a ,ifb) will be given by the 
condition 



J m (K a )d-(K h )J m+1 {K h ) - J m (^ b )d_(X a )J m+1 (^ a ) = 0. 



(24) 



Case ii) e > \h\: the two quantum-well states involved in the dot solution belong to different branches, e + and £_, 
with real values of K, noted K a and K^: 



* m { P M = c, { d+{Ka)Jm+l{KaP) s^) +c * { 



j m {K hP y m * 
d_(K h )j m+1 (K hP y( m+1 ^ l ' 



The boundary condition leads to 

J m {K & )d^{K h )J m+1 {K h ) - J m {K h )d+{K & )J m+1 {K & ) = 0. 



(25) 



(26) 



Case iii) —\h\ < e < \h\: one imaginary value of K, K a = in & , and a real value are involved. The energy associated 
to K a is 



and the coefficient for the wave function 



d±(K a ) = H) £± / a+ 5, = ~iS±(n a ) 
a' n K a + y«| 

With J m (iK P ) = i m I m (K P ), the quantum-dot wave function is then written as 



® m (p,<p) = i m c a 



i m ( KaP y m v 
s ± ( Ka )i m+1 ( KliP y ( - m + 1 ^ 



Cb 



j m (K hP y m ^ 
d^(K h )j m+1 (K hP y ( - m + 1 ^ 



The boundary condition leads to 

I m (K a )d_(K h )J m+1 (K h ) - J m (K h )S ± {K a )I m+1 (K a ) = 0. 

I 



(27) 
(28) 

(29) 
(30) 



Equations ( 24 ) , ( 26 1 , and ( 30 1 express a root-finding problem, which we solve numerically. We find a family 
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FIG. 2. The discrete eigenenergies of a quantum dot with radius R = 275 A and height L — 100 A for n = 1, shown as a 
function of the spin-orbit coupling strengths a and 7, without magnetic field (left side) and with B — 2.5 T (right side). States 
with \j z \ = 1/2, 3/2, 5/2, 7/5, and 9/5 are represented by different black, blue, green, red, and orange lines, respectively. In the 
central panels we keep 7 = and vary a from zero up to at r = 0.571 eVA, reported in [16]. Conversely, in the following curves 
(outside panels) a is fixed at a T , and 7 increases from zero to its final value of 7 r = 571.8eVA 3 . In-between the panels, the 
values of the quantum numbers l z associated to the nearby states at zero SOC are indicated. 



of solutions for each value of m that correspond to the 
discretized energies of the quantum dot. Moreover, the 
solutions give access to the wave numbers {K a , K^] such 
that we can determine the coefficients {c a ,Cb} from the 
boundary condition and the normalization of the in-plane 
wave functions (23), (25), and (29 1. All of these solutions 
carry a well-defined value of j z = m + 1/2, and in the 
absence of a magnetic field, the j z and —j z solutions are 
degenerate. 

The results for the energy levels are presented in Fig. [2] 
as a function of the SOC coupling strength. The states of 
different \j z \ are shown with different colors. To show the 
effect of the spin-orbit coupling, we start from the case 
of vanishing SOC in the center of the figure and increase 
the SOC strength up to the predicted values a r and 7 r 
corresponding to the left and right edge of the figure. 

Without SOC (a = 7 = 0, inner edges of the plot), 
the electronic states can be characterized by the or- 
bital angular momentum l z along the z-axis and the spin 
s = ±1/2, in addition to the total angular momentum 
j z = l z + s. The values of \l z \ corresponding to the 
states are indicated in the center of the figure. With- 
out magnetic field (left side), the states characterized by 
(l z ,s) = (±|Z a |,±l/2) are degenerate. In the presence of 
a magnetic field (right side), the Zeeman energy splits 
the levels corresponding different spin orientations. In 
the presence of SOC, the orbital angular momentum and 
the spin get mixed, l z and s cease to be good quantum 
numbers, and only the total angular momentum quantum 
number j z , shown by the different colors in Fig. [^charac- 
terizes the states. It can be seen states corresponding to 



the same \l z \ at zero SOC are split by the SOC according 
to the different values of \j z \. 

In order to discriminate the effects of the different SOC 
terms, we increase the SOC in two steps. We first con- 
sider the usual Rashba-like problem by setting 7 = and 
varying the linear coupling strength a from zero up to 
a r = 0.571 eVA. This situation is depicted in the inner 
part of Fig. [2] where the left side corresponds to the case 
of zero magnetic field and the right side to B = 2.5 T. 
The ensuing step is to fix a at a T and raise the value of 
7 from zero to 7 r = 571.8 eVA . The result is matched 
with the previous one and traced by the adjoining curves 
in the outer panels of the figure. It must be noted that 7 
determines not only the cubic-in-fc SOC coupling, but it 
also enters in the linear-in-fc coupling [cf. Eq. (15)]. Con- 



sequently, at the end of each curve we find the energy of 
the quantum dot for the corresponding a T and 7 r . 

Note the significant effect of 7 on the eigenenergies, 
that leads to much stronger energy changes than a alone. 
It brings, for example, the lowest pair of levels with 
j z = ±1/2 (lowest black curves) down to energies that 
are below Ef. Moreover, level crossings occur as a func- 
tion of 7, changing the order of the states in energy 
with respect to the case of vanishing SOC. This hap- 
pens mainly for the lowest energy states of a given \j z \ 
that are pulled down by the SOC below the higher en- 
ergy states with lower values of \j z \. We remark that 
the full range of eigenenergies that we consider has not 
been explored in previous studies, and that we explicitly 
include allowed energy values that lie within the gap of 
the two-dimensional dispersion relation of Fig. [T] that 



FIG. 3. Spin textures in cylindrical quantum dots with L — 100 A and R = 275 A. Left and right panels show results for the 
lowest and the second lowest states with \j z \ = 1/2. The arrows and colors indicate the spin orientation as a function of the 
position in the xj/-plane. Below the disks, the same data are shown for a linear cut through the center of the sample. 



is — \h\ < E\ < \h\. It can also be observed that the 
Zeeman splitting shrinks as the SOC increases, while the 
spin mixture brought about by the latter increases ac- 
cordingly. This indicates that the effective g-factor in 
quantum dots is affected by the SOC and depends on 
the geometry. 



B. Spin texture of the eigenstates 

We now investigate the properties of the quantum dot 
eigenstates. The spinor states of Eqs. (23 1, (25), and (29) 



contain the information on the spin texture of the dot 
states. Without SOC and in the presence of a magnetic 
field, even a very weak one, the states are spin polarized, 
and the spin texture of the one-electron states is uniform 
throughout the dot. The appearance of a non-trivial spin 
texture is therefore a signature of the SOC, and can be 
seen as the degree of mixing of the two spin components 
in an eigenspinor. To obtain the spin texture correspond- 
ing to a state, we compute the expectation value of the 
spin operator 



(<x)(r) = * T (r)erW(r) 



(31) 



for each spatial point r inside the quantum dot. Because 
of the separability of the wave functions (10), the spin 



orientation is independent of the longitudinal coordinate 
z. Moreover, the rotational symmetry of the dots around 
the z-axis imposes that the resulting spin orientations 
present the same symmetry. Therefore, their projection 
on the (^-direction vanishes, such that the local spin di- 
rection 



(er) (r) = f cos[/3(r)] + z sm[(3(r) 



(32) 



has only radial and z-components. The angle of the local 
spin orientation with respect to the xy-plane /3 depends 
only on the radial coordinate r. We construct the full 
eigenstate solution ^ with energy e as (23), (25), and 



( 29 ), depending on the value of e, with the corresponding 



j f a (e) and K\>{e) obtained from a numerical solution of 
the quantization conditions ( 24 ) , ( 26 1 , or ( 30 1 . 



In Fig. [3j we present two examples of spin texture in 
cylindrical quantum dots of length L = 100 A and radius 
R = 275 A, in the presence of the full linear and the cubic 
SOC terms with the coupling strengths a r and 7 r pre- 
dicted in Ref. [15] . The left panels show the dependence 
of the spin orientation on the position in the ccy-plane 
for the lowest energy states that have \j z \ = 1/2. This 
spin texture corresponds to one of the two sublevels in 
the lowest Zeeman doublet, shown in Fig. [2] by the two 
lowest black lines. The other of the sublevels, that are 
degenerate at B = 0, has spin orientations with the sign 
of the z-component reversed. The right panels show the 
spin texture for the next higher levels that are character- 
ized by \j z \ = 1/2, corresponding to the second pair of 
levels (black lines starting at \l z \ = 1 in Fig. [2]). 



C. Effective g-factor in quantum dots 

The effective g-factor is experimentally accessible, and 
it is thus a widely studied quantity. An example are the 
measurements of Ref. [I] , where the effective g-factor has 
been observed to depend on the dot size with absolute 
values that are reduced as compared to the bulk effective 
g-factor g* Ri —14.7 (value from Ref. |26j). In the exper- 
iment, the effective g-factor is extracted from the linear 
term of the magnetic-field induced energy splitting 



AE = \g eS n B B\ 



(33) 



of two states that are characterized by the same \ j z \ and 
degenerate in the absence of a magnetic field. According 
to this definition, each quantum-dot state has its own ef- 
fective g-factor, and we will focus on the effective g-factor 
of the ground state which is often the most relevant one. 
To calculate the effective g-factor we can use different 
approaches. The most direct way is to set the magnetic 
field strength to a small finite value, e.g. B = 0.1 T, and 
to calculate the difference between the two lowest dot en- 
ergies, using the procedure of Sec. |IV A| Alternatively, in 
order to avoid the finite value of the magnetic field, we 



can express the effective g-factor as 



5cff 



1 dAE 

[i B dB 



= 9 



de 
dh 



(34) 



in terms of the sensitivity de/dh of the quantized dot 
energy levels with respect to the magnetic field, at h — 



0. To determine this derivative, we proceed as in the 
case of Rashba SOC treated in Ref. 18J, and derive the 



quantization conditions ( 24 1 and ( 26 ) for negative and 
positive in-plane energy e, respectively. The resulting 
expression for the effective g-factor is 



, sga(e)u(K a ) + u(K b ) 



X - 



u(K a )u(K h ) 

Jm{K a )J m+ i{K b ) 



' C(Jfa, Kb) [2K a + sgn(e) U '(if a )]^ 1 + sgn(e)C(K b , [2K h - u'^b)]" 1 ' 
where we have defined the functions 

C(«a, K h ) = J m (K h )J' m+1 {K a ) + S gn(e)J' m {K & )J m+1 {K h ) (36) 



and u{K) = a'K-j'K 3 . We denote by J' m {K) and u'{K) 
the derivatives of the functions J„ 
K. 



and u with respect to 



The expression of Eq. ( 35 1 is a generalization of the re- 
sult of Ref. [TS] , and it reduces to the result given in Eq. 
(13) of that paper in the case 7 = of vanishing cubic- 
in-k SOC. In order to compute the effective g-factor us- 
ing the analytic expression ( 35 1 , we first determine the 



eigenenergies and the corresponding pair of wave- vectors 
.Ka and by solving numerically the quantization con- 
dition of Eqs. (24) and (261 at h = 0, and then evaluate 
(351 using the obtained values. In Fig. [4] we present our 
results for different dot dimensions with lengths rang- 
ing from 50 A to 200 A and radii from 150 A to 500 A. 
We have checked that a direct numerical evaluation of 
the level splitting from numerically calculated energies at 
small values of magnetic field B yields the same results as 
Eq. ( 35 ) . In the figure, the numerical data for g e s (black 
dots) is plotted as a function of the inverse effective di- 
mensionless linear in-plane spin-orbit coupling a 1 ^ 1 (see 
Eq. (15)). The data corresponding to different dot sizes 
approximately collapses on a single curve. While a plot 
as a function of a' shows the same data collapse, the 
presentation of Fig. [4] allows for greater clarity in the 
comparison with experiment. Such a single-parameter 
scaling shows that the dependence of the ground state 
effective g-factor g e ff on L and R is, at least within the 
range of explored sizes, to a good approximation given 
by a function of a' . Thus, the main mechanism giving 
rise to a size-dependence of the effective ground-state g- 
factor is the L-dependent renormalization of the effective 
linear coupling strength a' by the cubic SOC 7, and its 
scaling with R. 

For a fixed value of L, the renormalized linear-in-A: 
coupling strength a' is proportional to l/u E R. Since 
u E oc R~ 2 , we have a' cx R such that the effective linear 
coupling decreases with decreasing R. It can be seen in 
Fig. [4] that the value of g c s increases (in absolute value) 



to 
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0.75 
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FIG. 4. The calculated effective g-factors (full black circles) 
for cylindrical quantum dots of different length L and radius 
R, plotted versus 1/a' defined in (151, with the SOC param- 
eters from Ref. [16]. The data points are for dots with radii 
R from 150 A to 500 A, and lengths values L in the range 
between 50 A and 200 A. The crosses represent experimental 
data from Ref. !4, , obtained with a magnetic field perpen- 
dicular to the symmetry axis of the quantum dot. The blue 
horizontal line indicates the bulk effective value q* ~ —14.7. 



towards the bulk effective g-factor g* (blue line) as R and 
thus the effective coupling a' decreases. An increase in 
R leads to a larger a', and according to Fig. [2] the Zee- 
man splitting of the levels decreases as the SOC increases. 
The consequence is that |<7 e ff| diminishes. Conversely, for 
a given radius R, the increase in L leads to a decrease of 
the effective linear-in-fc coupling a' , with the result of an 
approach of g c g to g* . 

The effective cubic-in-fc coupling 7' oc l/u E i? 3 oc 1/R 
increases when R decreases, and a competition between 
a' and 7' can be expected. However, the spectrum of the 
lowest energies is related to small values of K and mainly 
dominated by the linear SOC (see Fig.[T]), at least for the 
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not too small values of R that we consider. The scaling of 
the results with a' leads to the conclusion that the main 
effect of the cubic Dresselhaus coupling 7 is the renor- 
malization of the effective linear-in-fc coupling, and that 
the impact of the effective cubic-in-fc coupling strength 
7' seems to be of minor importance. However, the above 
arguments are relevant for the case under study of not 
too small R and low-energy dot states. More important 
effects of the cubic-in-fc coupling 7' can be expected for 
the g-factor of excited states and in dots with very small 
R. 

Similarly to the results presented in Ref. |26j . where a 
zincblende Hamiltonian with adjustable parameters such 
as the energy band gap magnitude was used, we find 
negative values for the ground state g-factor of the dot. 
However, in our case small positive values do occur for 
short pillbox-shaped dots. In general, and similarly to 
the theoretical results for Rashba SOC'-^as well as the ex- 
perimental values of Ref. [I] (crosses in Fig. [4]), our effec- 
tive g-factors are of reduced absolute value as compared 
to the bulk effective g-factor g*. While the qualitative 
behavior and size-dependence of our results are clearly 
consistent with the data of Ref. [I] , a direct quantitative 
comparison cannot be made since in the experiment the 
magnetic field direction is not aligned with the symme- 
try axis of the dots. Also, while the effective g-factor has 
been measured for very different values of L, only a small 
range of radii has been covered in Ref. [I]. 



D. Phonon-induced spin relaxation rate 

An important quantity characterizing the usefulness 
of quantum dots for spintronics applications is the spin 
lifetime, which is generally limited by interactions with 
acoustic phonons. We consider here the most relevant 
situation for possible applications, namely, an electron 
initially in the higher sublevel \i) of the lowest Zeeman 
doublet which relaxes to the lower sublevel |/) due to the 
emission of a phonon. The rate T of this process can be 
calculated using Fermi's Golden Rule 



with Si(Q) being the deformation potential given by 



Hi(Q) = ~oAiy/Q, 



(39) 



where So is a bulk-phonon constant. The quantity A\ — 
2 y h Cx contains the mass density g and the sample vol- 
ume V. The deformation potential has the same form as 
in the case of a zincblende structure. The term Ai(Q) ac- 
counts for the piezoelectric contribution and upon intro- 
ducing spherical coordinates (Q,9 p ,ip p ) for the phonon 
momentum, it reads 

Al ( Q ) = Ai qT7^ a ° cos0 p ( h33 ~ hx sin2 9 p) ' (40) 

where h x = /i 33 — 2/i 15 — /i 31 . In general, h^j are bulk 
phonon constants and Ao = Aire/n, where K is the dielec- 
tric constant and e the electronic charge. We emphasize 
that 8 P is the angle between Q and the z-axis (defined 
as the c-axis of the wurtzite structure). 

The potential of the transverse phonon mode is given 

by 



tft(Q) = A t (Q)e^' r 



(41) 



with 



A t (Q) = A t -^A sin6 p (h 15 + h x cos 2 9 p ) . (42) 

We emphasize that in wurtzite lattices the transverse 
piezoelectric potential has only one term, while in 
zincblende lattices it has two. The matrix element in 



Eq. (37) can be factorized, and the rate can be written 



2?r 



£ \Mx(Q)\ 2 \(fW lQ - r ^)\ 2 n(Q)S(AE - hu> 



2tt 



Q,A 

(43) 

where M\ = ~i(Q) + iAj(Q) and M t = A t (Q). We 
first note that the modulus of the momentum is fixed 
by the (5-function. Concerning the integral over the 
electronic coordinates, we remark that both the ini- 



El(/|CMQ)|i)| 2 KQ) + l]5(Ai;-^ A ), (37) tial and the final states denoted by $H${p,<f>,z) 



Q,A 

where Q is the phonon momentum. The label A = {1, t} 
refers to the longitudinal and the transverse modes, re- 
spectively, and n(Q) is the Bose-Einstein phonon distri- 
bution with energy huj\ — hc\Q where c\ is the sound 
velocity of the corresponding mode. The energy differ- 
ence between the two electronic states AE = E$ — Et 
determines via the (5-function the energy of the phonons 
involved in the relaxation process. The potential U\(Q) 
comprises both the deformation and the piezoelectric 
contributions^^ for wurtzite lattice structures. For the 
longitudinal mode, we have 



'2/L sm{rnrz/L) have the same z-dependent 
factor. Therefore, the integral corresponding to the ma- 
trix element in (43) can be further split into two parts 



by using cylindrical coordinates, leading to 
\{f\e^\i)\ 2 = \Z(e p )\ 2 \T(e p , Vp )\' 



(44) 



The integral over z can be performed analytically, yield- 
ing 

\ryta \i2 2(27rn) 4 (l - cosq z ) 



ql [(2vm)2 - ql 



(45) 



tfi(Q) = [H 1 (Q) + iA 1 (Q)]e 



(Q r 



(38) 



where the definition q z = QLcos(9 p ) has been used. The 
other integral T(0p,<p p ) reads 
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T(0 p ,pp) = j &pp[(ciJ mf {Klp) + c{j ms {Klp)) (<iJ mi (Kip) + ciJ mi (Kip)) (46) 

+ (cfd£j mf+1 (Kfp) + 44 J ™f+i( K lpj) («4, + i(if:p) + 4^4,+iWp)) 

x / dip exp[i(nii — mf)<f] exp [iQsin(9p cos(ip p — <p)pR] ■ 
Jo 

The integral over ip can be easily performed by applying the Jacobi- Anger relation 

oo 

exp [ix cos <p] — i l Ji (x) exp . (47) 



Upon replacing Eq. (47) in Eq. (46) and carrying out the integration over ip, terms for all values of I vanish except 
the one with I = rn/ — rrii. The integral then results in 



d(pe i{ mi -m f )<p e iQ sin e B cos^p-^pfl _ ^ e i{ mi -m t ){<p v -K/2) ( pR Q sin p ) , 



(48) 



As it can be seen in Eq. (48), the complex exponen- 



tial becomes a common factor in Eq. (46), and leads to 
|T(0 p , </? p )| 2 = f(0 p ), which is not surprising, since the 
cylindrical symmetry is not broken by the phonon poten- 
tial. 

In addition to the determination of the values of 
Kb, and AE from numerically solving the quantization 
condition p4l, (pel, and pojj as in Sec. IIV AI the cal- 



culation of the spin relaxation rate still involves an inte- 
gral over p, and a subsequent integration over 9 P (since 



neither 



T(8^tp p )\ 2 nor |M A (Q)| 2 depend on (p p , cf. 



Eqs. (39), (40) and (42)), that we perform numerically 
as well. 

In the numerical evaluation of the relaxation rate, the 
parameters we use are g = 5900 kg/m 3 , c\ — 4410 m/s, 
c t = 2130 m/s, k = 15.15, and So = 5.8 eV, all taken 
from Ref. [2~T| . As the piezoelectric bulk constants for 
InAs nanowires having wurtzite structure have not been 
obtained so far from micr oscopic ca lculations, we follow 
the standard prescriptior l 27 ! 28 * 31 ! 32 ! of estimating them 
from the cubic structure by the use of the relations 
his = h 3 i = (-1/V$)hu and h 33 = (2/y/3)h u . For 
hi^ we use the value of the zincblende structure case 
(3.5 • 10 s Vm" 1 [21]), assuming that at least the order of 
magnitude of that value should be correct for the wurtzite 
case. It can be seen from Eqs. (40), (42), and (43) that 



the relaxation rate due to the piezoelectric contributions 
is proportional to /i 2 4 and thus not extremely sensitive 
to its precise value. 

We present in Fig. [5] the results for the spin relax- 
ation rate in a cylindrical dot of length L = 100 A and 
radius R = 275 A, as a function of magnetic field, as- 
suming that the initial and final eigenstates in the re- 
laxation process are the two lowest energy states (i.e the 
first Zeeman-split sub-levels). The black curves show the 
zero-temperature contributions of the different phonon 
potentials separately. It can be seen that the transverse 
piezoelectric mode yields the dominant relaxation rate 
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FIG. 5. The calculated contributions of different acoustic- 
phonon potentials to the spin relaxation rate as a function 
of magnetic field for cylindrical InAs quantum dots of length 
L — 100 A and radius R = 275 A, with wurtzite lattice struc- 
ture. The black curves correspond to the zero-temperature 
relaxation rate yielded by the piezoelectric transverse (TA- 
Piezo; solid), the piezoelectric longitudinal (LA-Piezo; dot- 
ted), and the deformation (LA-Defo; dashed) potentials. The 
green line shows the spin relaxation rate due to TA-Piezo at 
finite temperature T = 10 K. 



for magnetic field strengths below 1.25 T. Moreover, the 
relaxation rate strongly increases with the magnetic field 
strength. This is due to the Zeeman splitting that makes 
phonons of higher energies relevant where the density of 
phonon states is increased. Quite long spin lifetimes, 
of the order of 10 ns, occur for magnetic field strengths 
around 1 T, and much longer lifetimes are obtained at 
weaker magnetic fields. 

The temperature dependence enters solely through the 
Bose- Einstein distribution in Eq. ( 37 1 , such that the in- 



crease of the spin relaxation rate with increasing tem- 
perature can be easily obtained. The result for the dom- 
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mating TA-piezo mechanism at 10 K is shown in Fig. [5] 
(green solid line). 

Though it is clearly dominated by the piezoelectric 
contribution, the order of magnitude of the deformation- 
coupling caused spin relaxation rates we find is consis- 
tent with the one of Ref. [27], where the singlet-triplet 
relaxation for an InAs nanowire-based quantum dot was 
calculated. However, in that work, only the deformation 
coupling was taken into account and assumed to domi- 
nate. The same assumption was made in Ref. |28) . where 
the electron spin relaxation in a similar quantum dot was 
calculated. In both references, the supposed dominance 
of the deformation over the piezoelectric potential was 
justified by the fact that they considered small semicon- 
ductor nanostructures. As explained in Ref. [33], there is 
a competition between the two components that depends 
on the size of the nanostructure. For instance, the leading 
role of the piezoelectric coupling for weak magnetic fields 
has been reported^ for quasi-one-dimensional "cigar- 
like" quantum dots in GaAs nanowires with zincblcnde 
structure. For GaAs quantum dots, a crossing between 
the deformation and piezoelectric-induced rate as a func- 
tion of the magnetic field was found in Ref. [35] - In our 
case this occurs as well, though for lower values of mag- 
netic field than those observed in Ref. [35]. For InSb 
nanowires, numerical calcul ations show that the defor- 
mation potential dominates! 34 * 35 ! That domination has 
been assumed to be present in general, for all narrow-gap 
semiconductors.^ In contrast, in recent measurements on 
an InAs nanowire-based quantum dot j 3 -^ the piezoelectric 
coupling was crucial for the determination of the phonon 
spectrum. 

We find that for InAs, which has a larger band gap 
than InSb, but smaller than GaAs, the spin relaxation 
rate is mainly driven by the (transverse) piezoelectric 
phonon potential for magnetic fields below 1.25 T. Be- 
yond this value, the deformation seems to overcome the 
piezoelectric contribution, but our theory does not allow 
us to treat stronger magnetic fields. 



V. CONCLUDING REMARKS 

Motivated by recent experiments on spin-orbit qubits 
in InAs nanowire-based quantum dotsP we have calcu- 
lated the electronic structure of quantum wells and cylin- 
drical quantum dots with wurtzite lattice structure, tak- 
ing into account the linear and cubic Dresselhaus spin- 
orbit couplings and a weak applied magnetic field. We 
found analytical solutions for the energy levels of both 
types of structures. For quantum dots, we worked along 
the lines of the solutions previously found in the pres- 
ence of the linear Rashba spin-orbit coupling! 17 ! 18 ! Our 
obtained solution allowed us to explore the spin texture 
and the effective g-factor of the energy eigenstates, and, 
furthermore, we calculated the phonon-induced spin re- 
laxation rate in the ground-state Zeeman doublet. 

The energy dispersion of quantum wells shows a strong 



influence of the cubic Dresselhaus term. In quasi- 
two-dimensional structures, the cubic term leads to a 
thickness-dependent renormalization of the linear term 
which is of the same order of magnitude as the bare 
term. This renormalization produces a level crossing of 
all the subbands (within the range of validity of pertur- 
bation theory) that is being reported here for the first 
time. Another level crossing, which is a consequence of 
the competition between the linear and the cubic terms, 
and which has been reported earlier in wurtzite-structure 
materials, is also obtained here. Although our analyti- 
cal results are valid for all wurtzite-lattice materials, for 
concreteness we presented full numerical results for InAs 
quantum wells only. 

The quantum-dot eigenstates have been obtained as 
linear superpositions of degenerate quantum-well states, 
and their associated energies calculated as functions of 
the Dresselhaus coupling constants with and without a 
relatively weak Zeeman term. Again, a strong influence 
of the cubic term is observed in the energy levels. The 
spin texture of the energy eigenstates is shown to be 
qualitatively modified by the presence of the Dresselhaus 
spin-orbit coupling. Further studies and applications of 
these states as qubits should take into account this sel- 
dom discussed feature of the states. Our analysis of the 
effective ground-state g-factor shows a remarkable scal- 
ing collapse when the data are plotted as a function of 
the effective linear in-plane spin-orbit coupling a', that 
contains a size-dependent renormalization from the cubic 
Dresselhaus coupling 7. The obtained results are consis- 
tent with existing experimental dataP even though the 
magnetic field orientation is not the same. The scaling 
of the data indicates that other size-dependent mecha- 
nisms are of minor importance. However, the cubic in- 
plane coupling 7' is expected to become more relevant in 
the case of very small radius R and/or for the g-factor of 
higher excited states. 

Finally, we have calculated the acoustic-phonon- 
induced spin relaxation rate between the lowest Zeeman 
sub-levels as a function of magnetic field. The different 
rates arising from the longitudinal deformation, longitu- 
dinal piezoelectric, and transverse piezoelectric contribu- 
tions for wurtzite structures have been calculated and 
compared. While our results for the spin relaxation rate 
due to the deformation mechanism are consistent with 
those of Ref. [57], we find that, in contrast to what is 
usually expected for small nanostructures the trans- 
verse piezoelectric phonon potential gives the dominant 
relaxation rate, at least for our case of cylindrical dots. 
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